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Abstract 

This work sets the statistical affine shape theory in the context of real normed divi- 
sion algebras. The general densities apply for every field: real, complex, quaternion, 
octonion, and for any noncentral and non-isotropic elliptical distribution; then the 
separated published works about real and complex shape distributions can be ob- 
tained as corollaries by a suitable selection of the field parameter and univariate 
integrals involving the generator elliptical function. As a particular case, the com- 
plex normal affine density is derived and applied in brain magnetic resonance scans 
of normal and schizophrenic patients. 

1 Introduction 

The literature of matrix-variate distributions (real, complex, quaternion, octo nion) tells us 
about a great effort for obtaining separ ately topi cs that were noticed recently (^Di'az-Garcia| 
(|2009t ). lDi'az-Garci'a and Gutierrez- Jaim cz (200g),|D laz-Garcia and Gutierrez- J aimcz (2010)) 
can be derived under a general approach. 

In fact, many models and techniques are explored first in the real case, and then their 
extensions to the complex case are proposed joint with all the necessary mathematical tools 
for their develo pmen t . Fro m the last 60 years we can citate hun dred s of examples o f these 
extens ions, see Herz ( 1955 ) and James (1964), Muirhead ( 1982 ) and Khatri (1965), DavisI 
( 1980l ) and Ratnaraiah et al\ (2005), amon g many others examp les. 



In the statistical shape theory context, Drvden and Mardia (^9^)_gives an im portant 



summary of diverse techniques in the real case. By other hand. lMicheas et al\ (|2006[ ) studied 
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some of the topics in iDrvden and Mardia I (Il998l) for the complex case. 

There a number of techniques in shape th eory, we focus in t his p aper in the afiine 
approach. Goodall and Mardia (1993) (see also Diaz-Garcia et al. ( 20031 )) proposed an al- 
ternative system shape coordinates termed configuration or af fine coordinates, randoml; 



indexed by a matrix multivariate gaussian distribution. Then, iCaro-Lopera et al. I (|200 



aly 

m 



extended this theory by replacing the normality assumption with a matrix multivariate el- 
liptical law. Those works were studied in the real field, so if we follow the tradition of the 
literature, we can expect an extension to the complex case, for example, by studying the 
new jacobians, integrals and computations. 

Instead of this, we propose an unified approach for the statistical theory of shape, by 
studying the real, complex, quaternion and octonion cases in a simultaneous way. As we 
shall see in section [2l these four cases are formally termed, real normed division algebras. 
However as usual, this type of generalisation have a price, in this case we need some concepts 
and notation from the abstract algebra. 

For the sake of completeness, the case of the octonions is considered, but is important to 
highlight that many of the results for the octonion field, only c an be conjectured, because , 
many theoretical problems about these numbers remain open, see lDrav and Manogue (I1999I) . 
In fact, the rele vance of the octonions for understanding the real world is not clear at present. 



see 



BaezI (I2OO2I ). 



Section [5] reviews some definitions and notation on real normed division algebras, also, 
some concepts and integral properties of Jack polynomials and generalised hypergeometric 
function are given; then in section [3] a Jacobian with respect to Lebesgue measure for real 
normed division algebras is obtained and the main results follows as a consequence; section 
m gives the affine shape distribution for several particular elliptical laws; and finally, section 
[5] shows an application from the literature of shape in complex case. 



2 Preliminary results 



A detailed discussion of r eal normed division algebras may be found in iBaezI (l2002h and 
Gross and Richard^ ( 1987 ], and of Jac k polynomials and hypergeo metric functions in lSawyeil 
( 1997 ). Gross and Richards ( 1987 ) and Koev and Edelman ( 2006[ ). For convenience, we shall 
introduce some notations, although in general we adhere to standard notations. 

For us a vector space shall always be a finite-dimensional module over the field of 
real numbers. An algebra ^ shall be a vector space that is equipped with a bilinear map 
m : ^ X ^ ^ ^ termed multiplication and a nonzero element 1 € termed the unit such 
that m(l,a) — m(a, 1) — 1. As usual, we abbreviate m{a,b) = ab as ah. We do not assume 
^ associative. Given an algebra, we shall freely think of real numbers as elements of this 
algebra via the map w 1— >■ wl. 

An algebra 5^ is a division algebra if given a, 5 G with ab — 0, then either a = or 
6 = 0. Equivalently, is a division algebra if the operation of left and right multiplications 
by any nonzero element are invertible. A normed division algebra is an algebra ^ that 
is also a normed vector space with ||a6|| = ||a||||6||. This implies that 3^ is a division algebra 
and that ||1|| = 1. 

There are exactly four real finite-dimensional normed division algebras: real numbers, 
comp l ex nu mbers, quaternions and octonions, these being denoted generically as 5^, see 
Baej ( 2002 ). All division alge bras h a ve a r eal dimension of 1, 2,4 or 8, respectively, whose 
dimension is denoted by (3, see lBaezi (j2Q02l Theorems 1, 2 and 3). 

Let £^ „ be the linear space of all n x m matrices of rank m < n over ^ with m distinct 
positive singular values, where ^ denotes a real finite- dimensional normed division algebra. 
Let be the set of all n x m matrices over ^. The dimension of g^"^"' over K is /3mn. 
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Let A £ g^"X'™^ then A* = A denotes the usual conjugate transpose. 

The set of matrices Hi € s^d^ that H5;Hi = is a manifold denoted „, 

termed Stiefel manifold (Hi is also known as semi- orthogonal {P — 1), s emi-unitary (/? = 2), 
semi-symplectic (/3 = 4) and semi- e.cept^onal type (/3 = 8) matrices. see lP^^ d'Mano.uj 
(Il999l )i. The dimension of V^,„ over is [(3mn — m{rn — l)/3/2 — m]. In particular, V/^ ^ 



with dimension over 5R, [m(m + l)/3/2 — m], is the maximal compact subgroup il^(m) of 
£^ „j and consist of all matrices H G g-™^™ such that H*H = Therefore, il^(m) is the 
real orthogonal group 0{m) (/3 = 1), the unitary group U{m) {l3 = 2), compact symplectic 
group Sp{m) {(3 = 4) or exceptional type matrices Oo{ra) (j3 — 8), for 5^ = 5R, S) or O, 
respectively. 

Denote by ef„ the real vector space of aU S G 5^™x" such that S = S*. Let be the 
cone of positive definite matrices S G g'™^™. then is an open subset of 6(^. Over SR, 
(3^ consist of symmetric matrices; over £, Hermitian matrices; over Jo, quaternionic Hermi- 
tian matrices (also termed self-dual matrices) and over D, octonionic Hermitian matrices. 
Generically, the elements of &^ are termed as Hermitian matrices, irrespective of the 
nature of ^. The dimension of (3^ over 3? is [r7i(m — l)f3 + 2]/2. 

Let 2)^ be the diagonal subgroup of >C^j_,„ consisting of all D G g'™xm^ _ diag(cii, . . . , dm)- 
For any matrix X G g'"xm^ denotes the matrix of differentials (dxij). Fi nally, we 
define the measure or volume element (c?X) when X G g^™xn^ j^/^ or „, see lDimitrirJ 
(|2002^. 

If X G g'"xni then (dX) (the Lebesgue measure in ^'"x™) denotes the exterior product 
of the /3mn functionally independent variables 

n m /3 

(c?X) = /\ /\ where dxij = y/^ dx,-*^'*. 

1=1 j=i fc=i 

If S G 6f„ (or S G T^(m)) then (dS) (the Lebes gue measure in or in T^(to)) denotes 
the exterior product of the TO(m+ l)/?/2 functionally independent variables (or denotes the 
exterior product of the m{m — \)j3/2 + n functionally independent variables, if sn G 3? for 
all i = 1, . . . , m) 

m P 

i<.j k—1 
m m (3 

f\dsu f\ f\ds[f, ifs,, G3fi. 

i—1 i<j k—1 

Generally the context establishes the conditions on the elements of S, that is, if Sij G 3?, 
G C, G i3 or G D. It is considered that 

{dS) = \\dsf^\ds,.M\dsf. 

i<j k—1 i—1 i<j k—1 

Note that, the Lebesgue measure (c?S) requires that S G that is, S must be a non 

singular Hermitian matrix (Hermitian definite positive matrix). 

If A G 2)^ then (dA) (the Legesgue measure in D^) denotes the exterior product of the 
/3m functionally independent variables 



(dS) 



(dA) ^ f\ f\dX 



i=l k=l 
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If Hi e V^,_„ then 

n m 

(HJdHi) = /\ f\ h*dh,. 

i—l 

where H — (H1IH2) = (hi, . . . , h,„|h„i+i, . . . , h„) G il^(m). It can be proved that this 
differential form does not depend on the choice of the matrix H2 When m — 1; vf „ defines 
the unit sphere in 5^". This is, of course, an {n — dimensional surface in 5"". When 
m — n and denoting Hi by H, (H*(iH) is termed the Haar measure on il^(m). 
The surface area or volume of the Stiefel manifold „ is 



Voir 



(H^dHi) = 



2^m^mn/3 /2 



(1) 



where Tf^Ja] denotes the multivariate Gamma function for the space S^, and is defined by 
T^[a] = f etr{-A}|Ap-('"-i)'3/2-i(rfA) 

m 



where etr(-) = exp(tr(-)) 
Gross and Richardd (|l987l ). 



denotes the determinant and Re(a) > (m — l)/3/2, see 



Let (B) be the Jack polynomials of B = B*, c orresponding to t he partition k = 
(ki , ■ . . /cm) of fc, fci > • • • > fcm > with X)i!Li = sce lSawveij (|l997n and lKoev and Edelman 
(|2006h . In addition. 



pF^{ai,...,ap;bi,...,bg;B) = X! X! 



K [h]K,...,[bq]K 



kl 



defines the hypergeometric function with one matrix argument on the space of hermitian 
matrices, where [a]f denotes the generalised Pochhammer symbol of weight k, defined as 

m 

[af^=Y[{a-{^- 1)13/2),, 

2=1 

where ^(a) > (m — 1)0/2 — km . and ( a),; = a{a + 1 ) ■ ■ ■ {a + z — 1), see Gross and Richard^ 
(|l987l ). iKoev and EdelmanI (|2006D and iDi'az-Garci'al (|2009[) . 

Lemma 2.1. //X G £^ 



then 



HiGV; 



(tr(XHi))2'=(dHi) = ^ 



V2;fe cf(xx*). 



(2) 



See Dfaz-Garcia and Gutierrez- Jaimez 

Now, we use the complexification = + i©^ of ©^. That is, ©^^"^ consists of all 
matrices X e (3^ff)™x™ of the form Z = X + iY, with X, Y e ©^. We refer to X = Re(Z) 
and Y = Im(Z) as the real and imaginary parts of Z, respectively. The generalised right 
half-plane ^ = '^IL + i&j in 6^'^ consists of all Z £ ©^'^ such that Re(Z) G «p^, see 
Gross and Richardsj (|l987l p. 801). 



The next result generalises one given in IXu and Fang I (ll989l). lTeng et al. I (Il989l) and 
Caro-Lopera et al. I (|2009.) for the real normed division algebras, see iDiaz-Garcial (|2009( ): 
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Lemma 2.2. Let Z e ^ and U G 6f„. Then, 

f /t(trXZ)|X|°-("'-i)^/^-iCf(XU)(rfX) = ff"^"^''| |Z|-°Cf(UZ-^)7, (3) 

7 = / /i(z)z™+'=-idw < oo, (4) 
for Re(a) > (to — l)/3/2 — fcm, k = (^i, • ■ • km) a partition of k, ki > ■ ■ ■ > km > with 

= ^■ 

3 AfRne shape distribution 

Start with the following definition extended from lGoodall and Mardia I (Il993h . 



Definition 3.1. Two figures X G n -^i ^ ^kn have the same configuration, or 
affine shape, if Xi = XE + Ins* , for some translation e G ^ and a E G £^ j^. 

The configuration coordinates are constructed in the two steps summarised in the ex- 
pression 

LX = Y = UE. (5) 

The matrix U G £^ contains configuration coordinates of X. Let Yi G £^ and 

Y2 G C^K,^, with q = N -K -l>l, such that Y (Y^ Y^)*. Define also U = (I | V*)*, 

then V = Y2Y^^ and E Yi where L G ^ is a He l mert sub-matrix. 

Consider the following extension of Caro-Lopera et al. ( 2009L Lemma 8) for real nornied 
division algebras. 

Lemma 3.1. Let {Y^/'^f = F G qJ?^,. H G \ifi{K), and let E = pi/^H such that Y = 
Upi/^H. Then 

(dY) = 2~-^'|F|(«-i)^/2(dV)(dF)(H*dH). 
Proof. Let E = F^/^H, where E G -C^^^ , H G 'd^iK) and F^/^ ^ q}^^. Therefore 

E*E = H*FH. 

Then from Diaz-Garcia and Gutierrez- Jaimed (|2009l Lemma 2.2) 



(d(E*E)) = |H*H|(-^-i)'^/2+i(dF) = |Hp-i)'^+2(dF) = (dF). 
Now, by iDiaz-Garcia and Gutierrez- Jaime j ( 20091 . Lemma 2.15) 



(dE) = 2-^|E*E|'3/2-i((i(E*E))(H*dH) 
= 2--^|H*FH|''/2-i(dF)(H*dH) 

= 2--^|F|''/2-i(dF)(H*dH) (6) 



By other hand, observe that 



V ) ^ - ( VE 



And b y differentiating and computing the exterior product, see lDfaz-Garcia and Gutierrez-Jaimez 
(|2009l . Lemma 2.1), we obtain 



(dY) = |E|^«(dV)(dE), 
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but observing that |E| = iF^/^fjj ^ |F|i/2, we get 

{dY) = |F|'^«/2(dV)(rfE). (7) 
By replacing (|6]) into ^ the desired result is obtained. □ □ 

Now, recall that X G £^ „ has a matrix multivariate elliptically contoured distribution 
for real nor med division algebras if its density, with respect to the Lebesgue measure, is 
given by (see Diaz-Garcia and Gutierrez- Jaimez ( 20091) '): 



/x(X) 



;h{ty [S-i(X-/x)*0-i(X -/.)]} 



|5]|/3n/2|0|/3m/2 

where fi £ £^ S g e The function ft, : -t- [0, oo) is termed the generator 
function, and it is such that u^^^~^h{v?)du < oo. 

Su ch a distribution is denot e d by X ^ £^^^(n, S(g)0 h), for real case see Fang and Zhand 
and iGupta. and Vargal (|l993h and iMicheas et all (|2006l ) for complex case. Observe 
that this class of matrix multivariate distributions includes gaussian, contaminated normal, 
Pearson type II and VI, Kotz, Jensen-Logistic, power exponential, Bessel, among other dis- 
tributions; whose distributions have tails that are weighted more or less, and/or distributions 
with greater or smaller degree of kurtosis than the gaussian distribution. 

Now, we have the mathematical and statistics tools for establishing the main density. 

Theorem 3.1. Let X ~ 5^ .^^^(/ix, <8) &,h). Then the affine shape density is given 
by 



^/3ifV2r^ [/?(iv - 1)/2] 



E 



— y 

l)/2 ^ 

[/9(^-l)/2]f 



1 



[I3K/2] |S|W2|U*S-iU|'3(w-i)/2 ^ tW [K{N - l)/2 + t] ^ r! 



E 



C'^(U*f2S-iU(U*S-iU)-i)7, 



where 



1 



m/2t 



(8) 
(9) 



and S — LSL*, /x = i/.tx and = S ^fi&fj,*. 
Proof. Define 

LX0^/2 ^ LZ = Y = UE, 

where {&^^^)^ = 0. Thus Y - £^_i^^(/20"i/^ S g) I, /i) where S = LSL*, = i/xx- 
Therefore the density of Y is given by 



1 



/i{tr[S-i(Y - /x0-i/2)(Y - fxe-'^^)*]}. 



Making the factorisation Y = UE = UF^^^fj and by Lemma TS. 11 then the joint density of 
U, F and H is 

|p|/3(g+l)/2-l 



2if|S|/3K/2 

where n = S^VQ^V* 



tr (FU*S-iU + fi) + tr (^-20-i/V*S~^UFi/2jjj 

X {U*dU){dF){dY), 
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Expanding h in series of power, the joint density of U, F and H becomes: 



' ' . ^ [tr (FU*E-iU + n)] [tr (-20-^/ VS-^UFI/^H 



2if|S|/3K/2 



t=0 



(H*dH)(dF)(dV). 



Now, using the Lenima [2TT] for to integrate with respect to H, and recalhng that (l/2)t4*/ {2t)l — 
1/tl and C^(aB) = a*C^(B), the marginal joint density of F and U is 

1 ^1 y i./,(2t) rtr (FU*s-iu + n)] 

X V (U*f2S-iUF) (dF)(dV), 

Assuming that can be expanding in series of power, then the joint density of F and 

U is 



|S|W2r^ [f3K/2] ^0 H ^ ^ ^ ^ 

X y (U*nS-iUF) (dF)(dV). 



Hence, the marginal density of U is 
;r^^V2 



(10) 



|S|^^/2r^ [/3i^/2] ^0 H L V ;j .1. 

X /" (FU*S-iU)) |F|'9(?+i)/2-iC^ (U*f2S-iUF) (dF). 

Jf>o 

From ([3]), the integral in (|10p is evaluated as 

/ (FU*S-iU)] |F|^(«+i)/2-iC^ (U*rJS-iUF) (dF) 

Jf>o 

_ [^(A^ - l)/2]^ [/3(A - l)/2] C^(U*OS-iU(U*S-iU)-i) 
~ r[K{N -l)/2 + t] |U*S-iU|'3(^-i)/2 ^' 

where 

and the required result follows. □ □ 

Note that the general density is indexed by a simple univariate integral involving the 
general derivative of generator function . Thes e kind of densities appear rare in matrix- variate 
distributions (see Caro-Lopera et al. ( 2009f )). However, they demand the computation of 
derivatives of any order, which is not a trivial fact; general formulae for the classical elliptical 
models (Kotz, Pearson, Bessel, Jensen-Logistic) are available too in the above mentioned 
reference. 

Finally, the central and isotropic affine shape densities are obtained. 
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Corollary 3.1. 7/ X f^_-^j^^(0, Sx ® then the central affine shape density is 

invariant under the elliptical distributions, moreover, its density is 



4 [/?(iV-l)/2] 



Proof. The proof follows by taking f = r = in © and noting that h^'^^'^^^z) = h^°\z) 
h{z), hence 

|s|W2r^ [/3if/2] r[/?x(7v- i)/2] 

Now, using iFang and Zhand (Il990l p. 59), 



-l)/2 

7- 



,= / hiz)z^-i--^y^-^dz^'^^^^^''-'^/'^ 







^pK(N-l)/2 ' 

the desired result is obtained. □ □ 



Corollary 3.2. // Y ^ £^_^^^{^l^, a'^Ipf-i 'Si ®,h), then the isotropic noncentral affi 
shape density is given by 

.^^W^ [/3(A^ - l)/2)] - 1 yirtwn )V 

r^[/3if/2]|lK + V*V|^(^-i)/2^<!r[if(A^-l)/2 + t]^^r!^ 

where 



Tie 



and fli — a ^ fi& ^ fi* 



Proof. The result follows easily, just recall that U = (Ia' | V*)* and take S = cr^j^ 
dH). Thus 

□ 

4 The Gaussian affine shape distribution 

In this section we study the Gaussian affine shape distribution, as corollary of the preceding 
results. 

Corollary 4.1. Let X ^ ^N-ixKif^x.j S &,h). Then the affine shape density is given 
by 

[I3{N - l)/2] etr{-l3{n - U*rJS-iU(U*S-iU)-i)/2} 
TT^^i/^r^j^ [I3K/2] |S|'3^/2|U*S-iU|'9(^-i)/2 

iFf (-/3(j/2; /3if/2; -^U*nS-iU(U*S-iU)- 72). (11) 
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Proof. From iDiaz-Garcia and Gutierrez- Jaimez (l2009l) . the gaussian case turns, 

1 



h{v) 



27r\ 



PK(N-l)/ 



— e^^{-Pv/2} 



and 



/5 



2t+r 



27r\ 



Therefore 



/l(2*+-)(z)z/3^(^-l)/2+t-ld^ 



2t+r 



£13? f2TT\ 



pK(N-i)/2 cxp{-/3z/2}dz 



r[/3ir(iV- l)/2 + t] / /^/3 



TT' 



Also, note that 



[trO]-(-|) = [tr-/30/2]'- 



and 



^ (u*ns"iu(u*S"iu)-i) = c^(/3u*rjs-iu(u*s-iu)-V2)- 



The fi nal expression for (jlip is obtained by applying the Kummer relations, see lDiaz-Garcial 

(Hooi). □ □ 

The isotropic case of this distribution is given by 

Corollary 4.2. Let X ^ '^a'-ix/cI/^Xi &,h). Then the affine shape density is given 

by 

[I3{N - l)/2] etr{-/3(f2 - U*rJS-iU(U*S-iU)-i)/2} 



7r^^9/2r^ [f3K/2] \Ik + V*V|'3(^-i)/2 
li^f (-/3g/2; I3K/2; -/3U*f2S-iU(U*S-iU)- 72) 
[^(A^ - l)/2] eti{~(3{a-^fifi* ~ g-2-QH.^^H.u(^*^)-i)^2} 
7r/3W2r^ [PK/2] \Ik + V*V\0iN-i)/2 

iFf (-/3(7/2; (3K/2; -/3a-2u>/x*U(U*U)- V2). 



(12) 



(13) 



It follows straightforwardly, recall that U — (Ik \ V*)* and replace S = a^I^^i in 
Thus 

|s|W2|u*s-iu|''(^-i)/2 ^ ^ v*v|^(^-i)/2. □ 
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5 Example 



In this last section we apply the complex normal shape model in a classical data, th e 
brain magnetic reson a nce sc ans of normal and schizophrenic patients (see Bookstein ( 19661 ). 



Drvden and Mardia I (ll998l) . among many others). The random samples consist of 14 scans 



of each group, which correspond to a near midsagittal two dimensional slices of MR. On each 
image, an expert locates the following 1 3 landmarks and registers their co ordinates. The 



selected points correspond to ([Bookstein I (|l966l ). lDrvden and Mardia I (|l998f )): 1. splenium 



2. genu, 3. top of corpus callosum, 4. top of head, 5. tentorium of cerebellum at dura, 6. top 
of cerebellum, 7. tip of fourth ventricle, 8. bottom of cerebellum, 9. to p of pons, 10. bottom 



of pon s, 11. optic chiasm, 12. frontal pole, 13. superior colliculus (see lDryden and Mardia 



(|1998l . figure 9, p. 12)). The aim of the experiment is to test the equality in mean shape of 
the two population after removing the non geometrical information of the scans. 

The preceding works have analysed the shape under Euclidian transformation by filtering 
out translation, scaling and rotation of scans, and concluding that both mean shapes are 
statistically different. However, the shape theory via similarity transformation works well 
when the objects are rigid and develop a"constant radial growth" in some sense, but in the 
case of the brain, which is a soft organ, it is prone to deformations, so in order to match 
the scans we need to filter out the shear, instead of the rotation. Then the affinc shape or 
configuration analysis is more appropriate than the usual Euclidian shape ( this discrepancy 



can be solved statistically by using the so termed modified BIG criteria, see I Yang and Yang 



(|2007l l and the references therein). 

Now, if the original landmark 13x1 vector X follows an isotropic complex normal model, 
then from coroUarv 14.21 we have that the corresponding affine shape density is given by 

[Pi^ - l)/2] etr{-/3cr-VAt*/2 + l3a~^V* nn*V(lJ*lJ)-^ /2} 
7r/3^f9/2r^ [I3K/2] |U*U|^(^-i)/2 

iFf(-/3(?/2;/?i^/2;-/3a-2U>/x*U(U*U)-V2). (14) 

where /3 = 2, K = 1, 7V = 3, q = 11. Then it is of interest the estimation of the scale 
parameter cr^ and the location parameter 

M = + *Ml,2, ^J■2^ + A*3,1 + «M3,2, M4,1 + iM4,2, M5,l + «M5,2, ^J■6.1 + «/^6,2, 

/^7,1 + ^J■8A + «A'8,2, + «M9,2, MlO.l + iMlO,2 , + «A'll,2, Ml2,l + *A^12.2)' 

Note that the above affine density is a polynomial of degree eleven, then the inference 
procedure is notoriously simple. 

Let £(/i, a^, h) be the log likelihood function of a given group-model. The maximisation 
of the likelihood function £{'jl,a'^ ,h), is obtained in this paper by using the Nelder-Mead 
Simplex Method, which is an unconstrained multivariable function using a derivative-free 
method; specifically, we apply the routine fminsearch implemented by the sofware MatLab. 
The shape densities are polynomials of scalar zonal polynomials, i.e. hypergeometric se- 



ries w hich terminates and this can be computed easily by the algorithms of lKoev and Edelman 



(120061) 



At this point the log likelihood can be computed, then we use fminsearch for the 
MLE's. The initial value for the algorithm is the sample mean of the normal variables 
Y ~ A/i2xi(Mj'''^Ii2) and the median of the variances. 

The maximum likelihood estimators for location parameters associated with the normal 
and schizophrenic groups under the complex Gaussian model, are the following: 

For the normal group: 

fl = (-1.5312 + iO.6468, -0.0084 -I- iO.4637, 1.4306 + iO.9779, 1.2459 - i2.1091, 0.3962 - i0.7200, 
-0.0475 - il.3257, -0.2321 - il.9199, -0.6937 - i0.2605, -0.6556 - il.0118, -1.3664 + i0.1697, 

-2.4834 -I- il.4590, 0.2667 -|- i0.0804)', 
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and cr^ = 0.0159 and for the schizophrenic group: 

/I = (-1.6898 + j0.0545, 0.0179 + iO.5385, 1.1107 + il.3683, 1.7652 - il.7264, 0.5794 - iO.6484, 
0.3610 - il.300, 0.3015 - a.9547, -0.6052 - iO.4795, -0.3432 - il.2216, -1.3368 - jO.2235, 

-2.8060 + i0.6043, 0.1087 - i0.0434)' 

and 'P = 0.0155. 

Finally, we can test equality in afhne shape between the two independent populations. 
In this experiment we have: two independent samples of 14 patients and 24 population 
shape parameters to estimate for each group. Namely, if is the likelihood, 

where /x^, represent the mean shape parameters of the normal and schizophrenic 
group, respectively, then we want to test: Hq : /x„ = /x^ vs li\ : /x„ ^ fi^. Then 
—2 log A — 2sup^^ logL(/x„, /x^) — 2sup^jj log L(/x„, /x^), and according to Wilk's theorem 
—2 log A ~ X24 under Hq. 

In this case we obtained that: 

-2 log A = 2(858.09) - 2(834.60) = 46.98, 
Since the p-value for the test is 

Pixli > 46.98) = 0.0034 

we have important evidence that the norm al and schizophrenic brain M R are different in 
affine shape. This conclusion is ratified by iDrvden and Mardia I (|l998h . for example, but 



using a "rigid" Euclidian match in in this case they obtained a p-value of 0.005. 

Some new models (Kotz, Pearson VII, Bessel, Jensen-Logistic) can be studied and con- 
trast them with the Gaussian one via BIG modified criteria, however new extensions of the 
so termed generalised Kummer relation s with Jack polynomials are required (a generalisa- 
tion to real normed division algebras of iHerz (tl955h ): this shall constitute part of a future 



work. 

Finally, observe that the real dimension of real normed division algebras can be exp ressed 
as pq tentia of 2, /3 = 2" for n — 0,1,2,3. On the other hand, as observed by iKabd 



the results obtained in this work can be extended to hypercomplex cases; that is, for 
complex, bicomplex, biquaternion and bioctonion (or sedenionic) algebras, which of course 
are not division algebras (except the complex algebra), but are Jordan algebras, and all their 
isomorphic algebras. Note, too, that hypercomplex algebras are obtained by replacing the 
real numbers with complex numbers in the construction of real normed division algebras. 
Thus, the results for hypercomplex alg ebras are obt ained by simply replacing /3 with 2/3 i n 
our results (we reiterate, as reported bv lKabe (1984)). Alternatively, following Kabd ( 1984 ). 



we can conclude that, our results are true for '2^-ions', n = 0, 1, 2, 3, 4, 5, emphasising that 
only for n = 0, 1, 2, 3 are the result algebras in fact real normed division algebras. 
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